Surfactant-induced migration of a spherical drop in Stokes flow 
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In Stokes flows, symmetry considerations dictate that a neutrally-buoyant spherical particle will 
not migrate laterally with respect to the local flow direction. We show that a loss of symmetry due 
to flow-induced surfactant redistribution leads to cross-stream drift of a spherical drop in Poiseuille 
flow. We derive analytical expressions for the migration velocity in the limit of small non-uniformities 
in the surfactant distribution, corresponding to weak-flow conditions or a high-viscosity drop. The 
analysis predicts that the direction of migration is always towards the flow centerline. 

PACS numbers: 47.15.G-, 47.55.Dk 
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A classic result in microhydrodynamics is that lateral 
migration of a neutrally-buoyant, non-deformable spher- 
ical particle is prohibited in the creeping- flow limit. The 
result arises from the linearity of the Stokes equation and 
boundary conditions, and the symmetry of the problem 
under flow-reversal^. However, such cross-stream drift 
may occur if the symmetry is lost, e.g. by particle de- 
formation in a shear gradient or in the presence of a 
wal£. A small-deformation solution for a drop with a 
clean (surfactant-free) interface predicts lateral motion 
in unbounded Poiseuille flow, with the direction of mo- 
tion depending on the ratio of drop and suspending fluid 
viscosities^. In this note, we investigate the possibility 
of cross-stream migration of a non-deforming spherical 
drop, induced by asymmetric interfacial tension result- 
ing from redistribution of surfactant. 

Let us consider a drop with radius a and viscosity A77, 
embedded in an ambient fluid with viscosity 77. A surfac- 
tant, insoluble in the bulk phases, is adsorbed on the drop 
interface; diffusion of surfactant is neglected. The aver- 
age surfactant concentration is r eq , the corresponding 
interfacial tension is cr eq , and cto is the interfacial tension 
in the absence of surfactant. The drop is placed in a plane 
Poiseuille flow v°° = Uq ay 2 x , at position (xo, yo, Zo). In 
a coordinate system centered on, and translating with, 
the drop, the flow becomes 



Marangoni number 



= U, 
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ay 2 ) x - 
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where a is a measure of the curvature of the flow profile, 
7 = 2ayo is the local shear rate, and the migration ve- 
locity U m is the difference between the velocities of the 
drop and the undisturbed flow at the drop center. In 
creeping flows, a drop remains spherical provided that 
the capillary number is small 



Ca=^«l. 



(2) 



Convection by the surrounding flow creates non- 
uniformities in the surfactant distribution. The magni- 
tude of this redistribution is determined by the inverse 



Ma = — — ; Act = cto — c, 
Act 



cq j 
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which is the ratio of viscous stresses to the characteris- 
tic Marangoni stresses (gradients in surface tension) . For 
a dilute surfactant monolayer, a linear "perfect gas" in- 
terfacial equation of state relates surface tension and lo- 
cal surfactant concentration: cr(T) = ct(1) — Ma(T — 1). 
Henceforward, all quantities are normalized using a, 77, 
U , and r eq . 

By the linearity of the Stokes equations, the drop mi- 
gration arises as a superposition of that due to the flow 
about a clean drop, and that due to the flow driven by 
the Marangoni stresses 5 . The clean drop velocity field 
gives rise to slip, but no lateral migration 6 : 



U m = a- 



X 



3A + 2 



UL(T). 



(4) 



To determine the surfactant contribution U^, we 
solve for the velocity fields corresponding to the 
surfactant-covered drop using the formalism developed 
by Blawzdziewicz et al. 5 . In the original work, the 
method was applied to the dynamics of a stationary 
surfactant-covered drop in a linear flow. We have gener- 
alized the approach to treat a migrating drop in higher- 
order flows. Owing to the spherical symmetry of the 
problem, all quantities are represented in terms of spheri- 
cal harmonics ( Appendix [A"|) . Accordingly, the surfactant 
concentration T is expanded in scalar harmonics 



00 j 



j=l m=-j 
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The velocity field is expanded in a set of fundamental 
solutions of the Stokes equations 5,7 



(6) 



V ou t — c jmq [ U jm<j( r ) U jm(7( r )] + c jmq^j mq (r) , 
Vin = c j'ni(3U~ l m(2 (r) . 



Summation over repeated indices is implied. The ve- 
locity coefficients are decomposed into clean-drop and 
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surfactant contributions: Cj mq = cS + c| (T). The 
functions are vector solid harmonics related to the 

harmonics in the Lamb solution. The function 11^2 is 
radial, while and u^ ml are tangential to a sphere; 

u jmi i s surface-solenoidal (V s • u^ ral = 0). The far-field 
(imposed) flow is specified by v°° = c^ g u+ ; coeffi- 
cients for the Poiseuille flow {1} are listed in Appendix IB"1 
The velocity fields defined in ^ are naturally continuous 

across the interface because the basis fields u^L,„ reduce 

j in( i 

to the corresponding vector spherical harmonics yj mq at 
r = l. The velocity coefficients Cj mq are determined from 
the stress balance equations listed in Appendix O A 
fixed spherical shape limits the surface velocity to rigid 
motions. Thus, for j >1, Cj m 2 =0, and the other two co- 
efficients Cj mq are determined from the tangential stress 
balances alone, as use of the normal stress balance in 
conjunction with shape specification over-constrains the 
problem. Both the tangential and normal stress bal- 
ances are used to determine the c\ mq , which correspond 
to translational and rotational components of the flow. 
We obtain for the surfactant-driven flow 



c jm0 



V2 
(3A + 2) 3 



Mag jr 



(1 - Sjj) y/Jjj+T) 
A + 1 (2j + 1) 



Magj 



(7) 



0, 



'1/ 



(3A + 2) 3 



-Mag jm , 



where Ski is the Kronecker delta. The solution for j > 1 is 
identical to that of Blawzdziewicz et aL - . The clean-drop 
flow is 



u jm0 



3A + 2 
A + l 



(2A + 3) C ~ + V2(A-l)c~ 2 

3 r^y 



9r°° ~ 



2j + l 



(8) 



Jml 2 + j+A(i-l) jml ' 



u jm2 



3A 



2(A- l)c 



jmO 



(A + 4)c 



jm2 



The drop migration velocity is the difference between 
the volume-averaged velocity of the drop and the undis- 
turbed velocity at the drop center: 



U mig = -v°°(0) + J-^ v 



(r)dr . 



(9) 



Only the u+ m2 contribute to the integral, and only the 
yi m o and yi m 2 contribute to the first term. Hence, we 
obtain 



7 [-(Cll2 - Cl_l 2 )x 

-i(cn2 + ci_i 2 )y + V2c W2 z] 



(10) 



The above expression, in conjunction with ([8]) and (|B2[) . 
yields the slip term in (HJ) . Using ([7]) instead of ([5]) yields 
the surfactant-induced migration 



U s = 



Ma 



/67r(3A+2; 



[-foil - 9i-i)x 



-ifo 



ii + .gi-i)y + V2gwz] 



(11) 



The surfactant distribution coefficients gi m are deter- 
mined from the evolution equation^ 



dgjm 



Q t "J"' 1 L""J'"J2'"2 1 -*J"tJ2'"2J aj2'"2 (12) 

+ Ma[W(j)gj m + Qjmj 1 m 1 j2m29j2m 2 9jim 1 ] ■ 

The terms in this equation are defined in Appendix [D] 
The first three terms involve the clean-drop flow: Cj m de- 
scribes convection of the uniform part of the surfactant 
distribution by the imposed flow, Q describes rotation by 
the linear shear component of the flow, and A describes 
convection of the non-uniform part of the surfactant dis- 
tribution by the non-rotational components of the flow. 
The terms proportional to Ma pertain to flows driven by 
Marangoni stresses. The linear term describes relaxation 
towards the equilibrium uniform surfactant distribution, 
while the quadratic term describes convection of surfac- 
tant by the surfactant-induced flow. 

For arbitrary distortions of the surfactant concentra- 
tion, equation (|12p must be integrated numerically to 
determine U m . Here we consider steady-state solutions 
for the surfactant distribution, which implies that the 
surfactant evolution occurs on a faster time scale than 
the drop migration, i.e., r]a/Aa -C a/Um. Furthermore, 
we limit ourselves to small disturbances in the surfactant 
concentration, admitting an analytical solution for the 
migration velocity in the form of a regular perturbation 
expansion 



u -u (0) + u (1) 



(13) 



The choice of small parameter depends on the flow 
regime of interest. In weak flows, surfactant relax- 
ation towards the equilibrium distribution is fast, the 
gradients in surface tension are small, and Ma^ 1 is 
a natural choice. If the drop is far away from the 
centerline, yo 3> , the local flow is dominated by the 
shear 7 = 2ayo , and surfactant inhomogeneities are 
limited by drop rotation. In this case, the ratio of 
rotational and extensional time scales^ is ~A _1 , so this 
is an appropriate small parameter. Next, we examine 
these two regimes in more detail. 

Weak flow / Nearly-incompressible surfactant: 

Ma^ 1 < 1 and A = o{Ma) 

If MoT 1 <C 1, following the discussion in Blawzdziewicz 
et a/£, we introduce a regular expansion for the sur- 



factant concentration: gj m = ^ 



fe=0 



Ma 



-k-l(k) 



At 
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leading order, the evolution equation (fT2|) becomes 
= C r 
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(o) . . 2j + l 



WOOffjm- Solving for ^ yields 



JO) 



— ( r°° -\/2r° 
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and 



2VW 



l) c jm0 



3c 



jm2 



(14) 
for j > 1 . 



J(j + 1, 

Inserting the expression in (fTTj) and ([!]) gives 

UL 0) = ^x. (15) 

The Marangoni stresses immobilize the interface, the sur- 
face flow is incompressible, and the drop behaves like a 
rigid sphere. 

At next order, the equation is 



= [0.j mhm% + Aj m j 2 m 2 ] gf 2 L 2 + W{j)g 



III 



o (0) (0) 

™ 'jmj 1 m 1 j 2 m 2 9 'j 1 m 1 9 'j 2 m 2 ' 



(16) 



and the g^ give rise to a cross-stream migration velocity 



U« = -Ma^cfyo 



3A+ 17 „ 
9(A + 4) Y 



(17) 



This migration is directed towards the centerline for all 
values of the viscosity ratio A. Both the quadratic a 
and shear 7 components of the imposed flow distort the 
surfactant distribution, and interact with the distortions 
created by the other. In addition, the slip velocity 
interacts with the shear component. The drop migrates 
in the direction of decreasing shear magnitude. 

Far from the centerline / High viscosity drops: 

A" 1 < 1 and Ma = 0(A) 



If A 1 -C 1, we expand g jm = J2T=o^ 



-k-i (*) 

y jm 



and 

introduce Ma = A -1 Ma. The leading order equation for 
the surfactant distribution is the A — > 00 limit of (|12p . 
Only the first convection term, linear Marangoni term, 
and the rigid rotation within the f2 term are 0(1). Thus, 



= C jm + fl jmj m9j° m + MaW{j)g 



7 (0) 
jm 



(18) 



where C(jm) = \~ 1 C(jm) and W(j) — X^Wd) as 
A — > 00. The surfactant dynamics reflect both rotation 
and Marangoni relaxation balancing convection by the 
imposed flow. Solving for yields 



o (0) 



2a 



=F2Ma + 3iay 



Inserting in pT|) and (H)) and keeping terms up to order 
A -1 leads to both slip 



a 
3 



2a 
9A 



~ 2 
a Ma 



9A 9a 2 y 2 + <±Ma 



(20) 



and cross-stream drift 



a 2 y Q Ma 



3A 



9a 2 y 2 



4Ma 2 



(21) 



As Ma ^00, A— > 00, or yo — > 0, rigid sphere behavior is 
recovered; cross-stream motion vanishes, and the slip is 
given by the Faxen result (fT5|) . While the cross-stream 
migration velocity is still always directed towards the 
centerline, it is now a non-monotonic function of the 

distance therefrom, with maxima at yo = ±|^^. 

In conclusion, we have shown that the presence of small 
amounts of surfactant can significantly affect drop mo- 
tions in quadratic flows. In contrast with a clean, de- 
formable dropi, for which the direction of migration de- 
pends on viscosity ratio, a surfactant-covered spherical 
drop always migrates towards the flow centerline. This 
result is in agreement with those found for other particles, 
such as capsules^ and vesicles^, whose interfaces are gov- 
erned by Marangoni-like stresses. The present analysis 
may be extended to treat general surfactant redistribu- 
tions, and the drop trajectory dyo/dt = U miV [yo,T(yo)] 
determined by numerical solution of the surfactant evo- 
lution equation. 
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APPENDIX A: HARMONICS AND VELOCITY 
FIELDS 

Scalar and vector spherical harmonics are defined as 



Y — 

1 jrn — 



2j+l (j-m)\ 
4-7T (j+m)! 



(-l)T™(cos9)e 



_ 1 

Yjmo = [j (j + 1)] 2 r VnYj m , 
Vjmi = -if x y jm0 , 



(Al) 



Yjm2 



= vY, 



jm 



where the P™ are the Legendre polynomials, and Vn is 
the angular part of the gradient operator. The velocity 
basis functions are 



1 jm0 



\r 1 (2-j+ jr- 



) YjmO 



+ ^r-i\j(j + l)]i (l-r- 2 )y jm2 , 



U jml — T 1 yjml i 

^2 = 1^(2-^(^)4(1- 
+ \r-i (j + (2 - j)r- 2 ) y 3 



) YjmO 



ml 7 



(A2a) 
(A2b) 

(A2c) 
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x jmO 



V- 1 (-(j + l) + (j+3)r 2 )y jm0 



- [j (j + 1)]* (1 - r 2 ) y jm2 , (A3a) 

u j m i = r j yjmi , (A3b) 
u/ m3 = 3' J '" 1 (3+j)( i f i ) i (l-r a )yimo 

+ l^ 1 (j + 3 - (j + l)r 2 ) y jm2 . (A3c) 

APPENDIX B: IMPOSED FLOW 

The unbounded plane Poiseuille flow fl} seen by a mi- 
grating particle is represented as 

v°° = c°° u + (Bl) 



with coefficients 

! 4 "'" 

c 3±30 — =c "\/ ins ; 



3±32 — --"V 35 
c 3±10 — ztza T^ V 7 i c 3±12 — 15:01 5 V 21 ' 

271 



c 2±20 = T71 



u 2±ll - ~ a 3 V 5 

c?° ± io = T («f - 2(^ T iC/m,,)) v1F> 



C2±22 = T7M 



_ / 2ir 
c 101 — 1 7\/ "3" 



(B2) 



— —9TI J— r°° — — 277 /£ 
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APPENDIX C: HYDRODYNAMIC TRACTIONS 
AND STRESS BALANCES 

Neglecting isotropic contributions, the hydrodynamic 
tractions on a sphere may be expanded in harmonics 

thyd = Tjmqyjmq, with 



jmq 



Eg' C< jmq' i^q'q ^q'q) + C jmq'T qlq , 



T jmq — ^J2q' C jmq'T^ q . 



(CI) 



where T± q = T±, T % = ±(2j + 1), T* = T 3 W 



U— r, and 



Iu = ±(i + 3)-3' T a = ± ( 2 J + 1 ) 
T±=T±=0. 

The tangential stress boundary conditions are 

^!So " ^jmo = MaVj{j + T)g jm , 

out _ in _ r, 
jml jml w 5 

where a linear equation of state relating surface tension 
and local surfactant concentration is assumed: cr(T) — 
cr(l) — Ma(r— 1). The normal stress balance for a spher- 
ical drop requires 



(C2) 



_out in 



'jm2 ~ 'jmS — -2Mag jm (C3) 
For j > 1, we set Cj m2 = and do not use the normal 
stress balance condition. 



APPENDIX D: SURFACTANT EVOLUTION 
EQUATION 

The surfactant conservation equation on a moving in- 
terface is 

dT 

— +V,-(v.r)+r(vn)V-n = ) (Dl) 

where v s is the velocity component tangential to the 
surface. In the case of a sphere, the mean curvature 
V • n = 2 . Representing all quantities as harmonics and 
decomposing the velocity field into clean-drop and sur- 
factant contributions yields (fT2"|) . The clean-drop terms 
are 

Cjm = + 1)] 2 c" m0 — 2c° m2 ) 

fW™ 2 = [j(j + l)f 2 c° Jimil C Q1 , (D2) 
A- - - \i(n 4- 111* r° C 00 - 2 r° C 22 



Summation over ji and mi is implied, and the c°„ co- 
efficients refer to equations ©. Using equations (7|) to 
define the vector by 



the surfactant terms are 



(D3) 



W(j) = \j(j + l)]*W (j)-2W 2 (j), 



e 



\j(j + l)}*Wo(h)C u »-2W 2 (h)C 



(D4) 

i22 



The Clebsch-Gordan coupling coefficients C arc 
C 00 = B 

^ JiO'i + 1) + J"0" + 1) - iaO'a + 1) 



n n j 
ooo 



c 01 = B 



L?'0' + i)ii(ii + i)] 3 

h h 3 - 1 




(D5) 



(s + l)(s-2j 2 )(s-2j 1 )(a-2i + l) 



-f22 



Jl J2 J 




where s = j + ji + j 2 



D 



(-1)' 



(2i + l)(2-7i + l) (2i 2 + l) 



4tt 



Jl J2 J 

mi m 2 — m 



and I ^ 32 J \ ^ en0 ^ eg fj^Q Wiener 37-svmbol 11 
mi m 2 m ' 
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